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Abstract 

Glassy matter, as subjected to high shear rates, exhibit shear thinning : i.e., the viscos- 
ity diminishes with increasing shear rate. Meanwhile one prominent difference between the 
transport in micropores and that in macroscale is the (relatively) larger roughness observed 
inside micropores. As the pore size decreases, the surface-to- volume ratio increases and there- 
fore, surface roughness will greatly affect the transport in micropores. By treating the glass 
as a shear-thinning matter and using the rate-dependent model together with the boundary 
perturbation method, we can analytically obtain the transport results up to the second order. 
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1 Introduction 

In recent years, considerable effort was geared towards understanding how glasses respond to 
shear [1-3]. Phenomena such as shear thinning (i.e., the higher the shear rate is, the smaller 
is the flow resistance) and 'rejuvenation' are common when shear flow is imposed. Unlike 
crystals, glasses age, meaning that their state depends on their history. When a glass falls out 
of equilibrium, it evolves over very long time scales [4-7]. 

Meanwhile most of the classical solutions of contact problems, starting from the Hertzian case, 
rely on the assumption of nominally smooth geometries, which is reasonable at large enough 
scales. However, real surfaces are rough at the micro- or even at the meso-scale [8-10]. The 
role of surface roughness has been extensively investigated, and opposite conclusions have been 
reached so far. For instance, friction can increase when two opposing surfaces are made smoother 
(this is the case of cold welding of highly polished metals). On the other hand, friction increases 
with roughness when interlocking effects among the asperities come into play. This apparent 
contradiction is due to the effects of length scales, which appear to be of crucial importance in 
this phenomenon. 

In this paper we shall consider the transport of glassy matter in micropores which have radius- 
or transverse-corrugations along the cross-section. The glassy matter will be treated as a shear- 
shinning material. To consider the transport of this kind of glass (shear-thinning fluids) in 
microscopic domain, we adopt the verified transition-rate dependent model [4-5] which was used 
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to study the annealing of glass. 

We noticed that Cagle and Eyring [5] tried a hyperbolic sine law between the shear (strain) rate 
: ^ and (large) shear stress : r (because the relaxation at the beginning was steeper than could 
be explained by the bimolecular law) and obtained the close agreement with experimental data. 
They can obtain the law of annealing of glass for explaining the too rapid annealing at the earliest 
time. This model has sound physical foundation from the thermal activation process [4] (a kind 
of (quantum) tunneling which relates to the matter rearranging by surmounting a potential 
energy barrier was proposed therein; cf., Fig. 1). With this model we can associate the (glassy) 
fluid with the momentum transfer between neighboring atomic clusters on the microscopic scale 
and reveals the atomic interaction in the relaxation of flow with (viscous) dissipation. 
The outline of this short paper is as follows. Section 2 describes the general mathematical and 
physical formulations of the framework. In this Section, explicit derivations for the glassy flow 
are introduced based on a microscopic model proposed before [4]. The boundary perturbation 
technique [8-9] will be implemented, too. In the 3rd. and 4th. Sections, relevant results and 
discussion are given therein. 

2 Mathematical and physical formulations 

By the beginning of this century the concept of activation entropy was included in the model, 
and it was considered that molecules go both in the forward direction (product state) and 
in the backward direction (reactant state). The development of statistical mechanics, and later 
quantum mechanics, led to the concept of the potential energy surface. This was a very important 
step in our modem understanding of atomic models of deformation. Eyring's contribution to this 
subject was the formal development of the transition state theory which provided the basis for 
deformation kinetics, as well as all other thermally activated processes, such as crystallisation, 
diffusion, polymerisation, etc. [4-5] 

The motion of atoms is represented in the configuration space; on the potential surface the stable 
molecules arc in the valleys, which arc connected by a pass that leads through the saddle point. 
An atom at the saddle point is in the transition (activated) state [4]. Under the action of an 
applied stress the forward velocity of a (plastic) flow unit is the net number of times it moves 
forward, multiplied by the distance it jumps. Eyring proposed a specific molecular model of 
the amorphous structure and a mechanism of flow. With reference to this idea, this mechanism 
results in a (shear) strain rate given by 
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where V/, = A2A3A,V^ = A2A3A1, Ai is the perpendicular distance between two neighboring layers 
of molecules sliding past each other, A is the average distance between equilibrium positions in 
the direction of motion, A2 is the distance between neighboring molecules in this same direction 
(which may or may not equal A), A3 is the molecule to molecule distance in the plane normal 
to the direction of motion, and r is the local applied stress, 5G is the activation energy, h is 
the Planck constant, ks is the Boltzmann constant, T is the temperature, Vh is the activation 
volume for the molecular event [4-5]. The flow of the glass is envisaged as the propagation of 
kinks in the molecules into available holes. In order for the motion of the kink to result in a 
plastic flow, it must be raised (energised) into the activated state and pass over the saddle point. 
Solving Eqn. (1) for the force or r, one obtains: 

r = -^smh^(|), (2) 
which in the limit of small {^/B) reduces to Newton's law for viscous flow. 

We start to consider a steady, fully developed flow of the glassy fluid in a wavy-rough microan- 
nulus of r2 (in mean-averaged outer radius) with the outer wall being a fixed wavy-rough surface 
: r = r2 + e sm{k9) and r\ (in mean-averaged inner radius) with the inner wall being a fixed 
wavy-rough surface : r = ri + esui(k9 + (3), where e is the amplitude of the (wavy) roughness, 
and the wave number : k = 2Tr/L {L is the wave length), /3 is the phase shift. The schematic is 
illustrated in Fig. 2. 

Firstly, this glassy fluid can be expressed as ^ = ^0 sinh(r/To), tq = 2kBT/Vh, where ^ is the 

shear rate, r is the shear stress, and 

/ _ o ^VhksT -6G 
^° = ^ = '^^^^Pfe^ 

is a function of temperature with the dimension of the shear rate (for small shear stress r ^ tq, 

tq/^o represents the viscosity of the material). In fact, the force balance gives the shear stress 

at a radius r as r = —{r dp/dz)/2. dp/dz is the pressure gradient along the flow (or tube-axis : 

z-axis) direction. 

Introducing the forcing parameter 11 = —(r2/2ro) dp/dz then we have ^ = ^0 sinh(nr/r2). As 7 = 
—du/dr {u is the velocity of the fluid flow in the longitudinal (z-)direction of the microannulus), 
after integration, we obtain 

■u = -I- ^^[coshn — cosh(— )], (3) 
n r2 

here, Ug is the velocity over the surfaces of the microannulus, which is determined by the bound- 
ary condition. We noticed that a general boundary condition for fluid flows over a solid surface 
was proposed (cf., e.g., [10-11]) as 

Sn = L^^l - i)-^/^ (4) 
sc 
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where 6u is the velocity jump over the solid surface, is a constant slip length and is 
the critical shear rate at which the slip length diverges. The value of is a function of the 
corrugation of interfacial energy. We remind the readers that this expression is based on the 
assumption of the shear rate over the solid surface being much smaller than the critical shear 
rate of (c- Cc represents the maximum shear rate the fluid can sustain beyond which there is 
no additional momentum transfer between the wall and fluid molecules [10-11]. How generic 
this behavior is and whether there exists a comparable scaling for glassy fluids remain open 
questions. 

With the boundary condition from [10-11], we shall derive the velocity field and volume flow 
rate along the wavy-rough microannuli below using the boundary perturbation technique (cf., 
e.g., [8-9]) and dimensionless analysis. We firstly select La = r2 — ri to be the characteristic 
length scale and set 

r' = r/La, Ro = r2/La, Ri = ri/La, e' = e/La. 

After this, for simplicity, we drop all the primes. It means, now, r, Ro, Ri, and e become 
dimensionless. The walls are prescribed as r = Ro + esm{k9), r = Ri + esm{kd + (5) and the 
presumed fully-developed flow is along the z-direction (microannulus-axis direction). Along the 
boundary, we have 

^ = ((iit /(in) I surface' 

where, n means the normal. Let u be expanded in e : u = uq + eui + e^U2 + • ■ ■, and on the 
boundary, we expand u(ro + edr, 6{= 9o)) into 

uir,9)\(^ro+€dr,eo) = u{ro,0) + e[drur{ro,9)] + e^[—urr{ro,e)] H = 

K., + ^cosh(g)|^+-^^(^^), (5) 

where the subscript means the partial differentiation (say, Ur = du/dr) and 

■^^siiplon surface ~ -^s^ii^ ~ T") ^ ]lon surface' (*^) 

sc 

nsu,, = i^^eo[sinhn(l - |o^)-i/2]. (7) 

Sc 

Now, on the outer wall (cf., e.g., [8-9]), 

■ du „ V (r — Ro — esm( k6)) 9/,/>M-ir i 

k 



5 



k 1 k"^ 

^uo,\R,cos{ke)] + e^[---^ 



1 A; 
7;U0rrr\Ro sin^ {kO) - cos{ke){ui,\R^ + uqJr^ sin(A;^))] + O(e^). (8) 



Considering ~ Ro,Ri ^ e case, we presume sinhll <C ^c/Co so that we can approximately 
replace [1 - (^o sinhn)/4]-V2 by [1 + sinhn/(2^c)]- With equations (5), (6) and (8), using 
the definition of ^, we can derive the velocity field up to the second order. The key point is to 
firstly obtain the slip velocity along the boundary or surface. 

After lengthy mathematical manipulations and using (1 — C/Cc)~^^^ 1 + ^/{2^c), 
uo = -^^[cosh(^) - coshn] + Usiip^, ui = sin(A;6') sinhll + UgHp^, 



we have 



■^Sc Sc -^sc Sc 



Uslipo + e '^slipi + e'^Uslip2 + O(e^) (9) 

where 

k n ■ n 

Uf = ui^ + uo^^ sm{k9) ^ cos{k9)uog = -^Co cosh(— r) sin(fc6'), 



Ro Ro 



and 



^isc = cos'^{k9)uo^ + -tto,,, sin^(A;e) = ^^oi-^ cos^(/c^) - sin^(A:e)] sinh(— r). 

Thus, dX r = Ro, up to the second order, 

Usiip = us = L% sinhn(l + ^) + eio sin(A:e)[sinhn + coshH (1 + Kq)] + e^L^,^{[ 

2 Kg 2 

„ ' sm\ke)-^ cos\ke)+-. sm\ke)] sinhn(l+ifo)+75^4^ cosh^ H sm\k9)}, (10) 

-Ko-t-s Ro -n-o -f^o 

where Kq = l + (^osinhn)/^c Prom the velocity fields (up to the second order), we can integrate 
them with respect to the cross-section to get the volume flow rate {Q, also up to the second 
order here). Q = J^" I^°+^s^lkeli3) ^(^' d)rdrd9 = Q smooth + eQpo + Qp2- 

Q = ^io{Ll{Rl - Rl) sinhn (1 - + ^[(^^2 _ ^^3i^n - ^RI sinh$- 

sinh(n^) + ^(coshn - cosh(n^)]} + e^-KUTl^^^iRo - |^)]+ 
eosinhn(l + f^)(-A;^ + )[1 - (|i)2] + f [(«.Hpo+ 
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^ coshn) + ioRoi- sinhn + + URi ^^H^§^) - ^ cosh(ng))]+ 

1 rr , „ „coshn, sinhn, „ . ^ . ^m-, 
7r,eo{[sinhn + H— — (1 + -^-^)][Ro - i?i(cos/3 + sin /?)]}+ 

4 -"o 



Here, 

usiip, = L°7o[sinhn(l - ^!^)-V2]. (12) 

3 Results and discussion 

In the following, as our interest is about the transport of glassy matter at very low temperature 
environment (H — ^ 0), thus we should take the asymptotic limit of Q. In fact, we have, for the 
smooth surface of the micropore, 

Qsmoothln^o = TrioiL'Ml - Rl) sinhn (1 - ^^)-i/2 + ^[^r^ _ coshn- 

^(Rl sinhn - i?,i2<,sinh(n^)) + ^(coshn - cosh(n^))]}|n^o = 0, (13) 

11 Kg 11 Ko 

and for the wavy-rough surface (up to the second order) Qp2 = as n — >^ 0. Now, the (referenced) 
shear-rate 

is a function of temperature and the activation energy. From equation (11), we can ob- 
serve that the transport of glassy matter in very-low temperature : T ^ is rather small 
(Qp2 oc $,0 (X rexp(— 1/T))! It is worth pointing out that the Eyring model requires the interac- 
tion between atoms in the direction perpendicular to the shearing direction for the momentum 
transfer [4-5]. This might explain why our result is orientation dependent (due to /3). 
Note that, based on the rate-state Eyring model [4] (of stress-biased thermal activation), struc- 
tural rearrangement is associated with a single energy barrier (height) E that is lowered or raised 
linearly by an applied stress a or r. If the transition rate is proportional to the shear strain rate 
(with a constant ratio : Co), we have 

a = E/V* + (kBT/V*) Hi/Couo), (15) 

where V* is a constant called the activation volume, ks is the Boltzmann constant, T is the 
temperature, and i^o is an attempt frequency [4-5,12]. Normally, the value of V* is associated 
with a typical volume required for a molecular shear rearrangement. Thus, the nonzero flow rate 
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(of the glass) as forcing is absent could be related to a barrier-overcoming or tunneling for shear- 
thinning matter along the wavy-roughness (geometric valley and peak served as atomic potential 
surfaces) in annular micropores when the wavy-roughness is present. Once the geometry-tuned 
potentials (energy) overcome this barrier, then the tunneling (spontaneous transport) inside 
wavy-rough annular micropores occurs. 

We also noticed that, as described in [4], mechanical loading lowers energy barriers, thus facil- 
itating progress over the barrier by random thermal fluctuations. The simplified Eyring model 
approximates the loading dependence of the barrier height as linear. The linear dependence will 
always correctly describe small changes in the barrier height, since it is simply the first term in 
the Taylor expansion of the barrier height as a function of load. It is thus appropriate when the 
barrier height changes only slightly before the system escapes the local energy minimum. This 
situation occurs at higher temperatures; for example, Newtonian flow is obtained in the Eyring 
model in the limit where the system experiences only small changes in the barrier height before 
thermally escaping the energy minimum. As the temperature decreases, larger changes in the 
barrier height occur before the system escapes the energy minimum (giving rise to, for example, 
non- Newtonian flow). In this regime, the linear dependence is not necessarily appropriate, and 
can lead to inaccurate modelling. This explains why we should adopt the hyperbolic sine law 
[4-5] to treat the glassy matter. 

Finally, wc present the calculated maximum velocity (unit : m/s) with respect to the temper- 
ature in Fig. 3. Geometric parameters : r2 = 100 nm, the activation volume : 0.2 nm^ and 
the roughness amplitude e = 0.05r2. We consider the effect of the activation energy : 10^^ and 
2 X 10^^ Joule. Around T ~ 0.25°K, the maximum velocity (of the glassy matter) either keeps 
decreasing as the temperature increases for larger activation energy or instead increases as the 
temperature increases for smaller activation energy! The latter observation might be related 
to the argues raised in [13] for annealing process of solid helium at similar low temperature 
environment if we treat the solid helium to be glassy at low temperature regime. 

4 Conclusion 

To conclude in brief, we obtain the transport of glassy matter inside annular wavy-rough mi- 
cropores under very low temperature environment. The flow rate is rather small (of the order 
of magnitude of the square of the small wavy-roughness amplitude) and is proportional to the 
(referenced) shear rate (which is strongly temperature as well as activation energy dependent), 
the slip length and phase shift of the wavy-roughness as illustrated above. Acknowledgement. 
The author stayed at the Chern Shiing-Shen Institute of Mathematics, Nankai University around 
the beginning of 2008-Jan. Thus the author should thank their hospitality for the first stage 
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Fig. 1. The structural contribution to the shear stress is : Shear thinning. Increasing 
strain causes a local energy minimum to flatten until it disappears (removing of energy barrier 
or quantum-like tunneling). 




Fig. 2. Schematic of an annular micropore. /3 is the phase shift between the outer 
and inner wavy-roughness, e is the amplitude of small wavy-roughness. 



10 



Activation voiume=0.2 nm 



Activation energy=10" Jouie 



Activation energy=2 x 1 Joule 



0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Temperature (°K) 



0.45 



0.5 



Fig. 3. Comparison of calculated (maximum) velocity (unit : m/s) using two activation 
energies 10"'^^ and 2 x 10~^^ Joule. Around T ~ 0.25°K, the monotonic trend of velocity 
bifurcates as the temperature increases. r2 = 100 nm and e = 0.05r2 here. 



